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program  is  the  command  file  described  above. 

The  remaining  three  programs  require  no  interaction  and  may  be  used  in  a  batch  mode.  Each 
of  these  produces  error  messages  whenever  unexpected  conditions  are  encountered.  Insofar  as 
possible,  these  messages  were  intended  to  be  self-explanatory.  If  more  details  are  required,  each 
program  also  contains  a  complete  list  of  error  messages  and  an  explanation  of  each. 

This  report  includes  descriptions  of  the  programs,  examples  of  the  required  inputs,  and  copies 
of  typical  program  outputs.  Complete  source  listings  are  also  provided  in  the  appendices. 
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The  METEOR  Software  Package  for  Analysis  of  Meteorological  Data 
I.  INTRODUCTION 

A  requirement  for  many  experiments  in  environmental  chemistry  is  that 
extensive  meteorological  records  must  be  maintained.  Frequently,  simultaneous 
collection  of  several  different  types  of  data  (temperature,  humidity,  pres¬ 
sure,  etc.)  at  regular  intervals  (minutes  to  hours)  over  rather  long, 
continuous  periods  (days  or  weeks)  is  necessary.  The  use  of  modern  data 
loggers  greatly  simplifies  the  acquisition  of  data  and  often  provides  for 
storage  and  transfer  on  computer-readable  media,  such  as  magnetic  tape.  Off¬ 
line  analysis  involves  conversion  of  the  data  to  physically  meaningful  units, 
calculation  of  derived  quantities,  and  presentation  of  the  results  in  formats 
which  are  convenient  for  the  user's  purposes.  Due  to  the  large  quantities  of 
data,  these  steps  are  usually  very  time-consuming. 

METEOR  is  a  FORTRAN  program  package  which  is  designed  to  alleviate  many 
of  these  difficulties.  It  provides  software  for  inputting  data  files,  search¬ 
ing  out  relevant  portions  of  these  files,  processing  data,  and  generating 
printed  or  plotted  output. 

The  goal  has  been  to  maximize  the  generality  of  the  program  while  mini¬ 
mizing  the  demands  on  the  operator.  The  former  goal  has  been  addressed  by  the 
liberal  use  of  subroutines  and  functions,  making  program  expansion  and  alter¬ 
ation  a  simple  matter  of  inserting  new  or  updated  program  modules. 

The  requirement  that  the  package  be  easy  to  use  has  led  to  the  develop¬ 
ment  of  an  auxiliary,  interactive  program  which  aids  in  the  creation  of  a 
command  file.  This  file  then  controls  the  execution  of  the  primary 
programs.  Currently,  there  are  three  primary  programs  (METSRT,  METCLC,  and 
METPLT)  and  one  auxiliary  program  (METINP)  in  the  METEOR  package.  The  rela¬ 
tionship  among  these  programs  is  illustrated  in  Figure  1. 

METSRT  (METeorological  data  SoRTing  program)  is  responsible  for  input  of 
raw  data,  selection  of  the  data  required  for  analysis,  preliminary  processing 
(conversion  of  units  and  scaling,  where  necessary),  testing  of  the  data  for 
various  error  conditions,  listing  of  selected  data  in  tabular  form,  and  cre¬ 
ation  of  a  file  containing  all  processed  data. 

METCLC  ( METeorological  data  CaLCulation  program)  then  operates  on  the 
output  of  METSRT,  calculating  the  values  of  various  derived  quantities  (total 
moisture  loading,  for  example).  These  results  are  added  to  the  file  of  pro¬ 
cessed  data  and  may  also  be  printed  out,  again  in  tabular  form. 

METPLT  ( METeorological  data  PLoTting  program),  the  last  of  the  primary 
programs,  reads  the  data  file  produced  by  the  previous  programs  and  generates 
selected  graphs. 

It  is  possible  to  string  these  programs  together  and  to  process  the  data, 
from  raw  data  input  to  finished  plots,  in  one  continuous  batch  job.  However, 
the  sequential  design  of  the  programs,  with  intermediate  outputs,  was  intended 
to  allow  easy  operator  intervention  in  the  event  that  bad  data  is  encountered. 
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Figure  1 

Relationship  of  the  programs  and  files 
in  the  METEOR  program  package. 
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The  auxiliary  program,  METINP,  (METeorological  parameter  INPut  program), 
is  interactive  and  is  intended  to  be  run  from  a  CRT  terminal.  This  program 
requests  information  reqarding  the  specific  functions  in  METSRT,  MET'CLC,  and 
METPLT  that  are  desired  and  constructs  a  control  file  in  the  appropriate 
format.  This  file  is  displayed  line-by-line  for  operator  verification  before 
being  written  onto  the  disk.  Errors  may  easily  be  corrected  at  this  point. 

The  METEOR  package  is  currently  running  on  a  DEC  System-10  and  the  1/0 
file  names  discussed  below  are  those  used  by  the  TOPS-IO  operating  system.  No 
major  problems  are  expected  to  arise  if  the  programs  are  transfered  to  another 
system  (in  fact,  earlier  versions  of  METSRT  and  METPLT  were  run  on  a  Texas 
Instruments  Advanced  Scientific  Computer).  However,  METINP,  because  it  is 
interactive,  will  not  operate  properly  in  a  batch  processing  environment. 

II.  METSRT  OPERATION 

For  reference  during  the  following  discussion,  a  complete  listing  of  the 
METSRT  source  code  is  given  in  Appendix  A. 

METSRT  may  be  logically  divided  into  input,  processing,  and  output  sec¬ 
tions.  The  input,  and,  to  some  extent,  the  processing  sections  of  the  program 
must  be  tailored  to  the  characteristics  of  the  data  source.  METSRT  was  origi¬ 
nally  written  specifically  for  use  with  a  Fluke  Model  2240B  data  logger, 
having  the  analog  data  output  format  shown  in  Table  1 .  A  digital  data  format 
(Table  1 )  is  also  available,  but  is  not  currently  in  use.  Other  formats  would 
necessitate  changes  to  the  search  and  input  routines  and  possibly  to  the  error 
flagging  subroutine.  We  are  presently  making  alterations  in  order  that  the 
output  of  a  newly  constructed  data  acquisition  system  may  be  processed  by 
METSRT. 

Four  I/O  files  (and  four  different  logical  devices)  are  involved: 

1 )  FOR04.DAT  (device  4)  contains  parameters  which  control 
program  execution. 

2)  FOR05.DAT  (device  5)  is  the  input  data  file. 

3)  FOR06.DAT  (device  6)  is  a  tabular  output  for  printing. 

4)  FOR07.DAT  (device  7)  contains  all  data  and  parameters 

and  is  intended  to  be  read  by  subsequent  programs. 

Initially,  METSRT  reads  the  control  file  (FOR04.DAT) ,  which  specifies  the 
dates  of  interest,  the  specific  types  of  data  which  are  to  be  processed,  the 
desired  format  for  printed  output,  and  the  units  for  both  input  and  output. 
Print  switches  may  be  set  to  select  channels  for  which  data  is  to  be  listed. 

In  addition,  other  parameters,  pertaining  to  METCLC  and  METPLT,  may  be  pres¬ 
ent.  These  parameters,  if  present,  are  ignored  by  METSRT.  An  example  of  the 
control  file  is  given  in  Table  2. 

A  subroutine  (SEARCH)  is  then  called  which  searches  the  data  file 
( FOR05.DAT)  for  the  first  data  set  and  reads  the  time  and  date  header.  Any 
data  set  which  is  dated  prior  to  the  specified  initial  time  is  rejected  and 
the  search  continues  until  one  of  the  following  conditions  is  met: 

1)  An  End  of  File  (EOF)  is  read. 
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Table  2 

parameter  file  (FOR04.DAT)  required  to  process 
a  portion  of  the  data  from  the  1981  cruise  of  the 
USNS  HAYES. 
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2)  A  data  set  within  the  specified  time  window  is  found. 

3)  A  data  set  having  a  date  later  than  the  last  desired 

time  is  encountered. 

Cases  1 )  and  3)  cause  appropriate  error  messages  to  be  printed  and  execu¬ 
tion  terminates.  It  is  assumed  that  time  monotonically  increases  between  data 
sets. 


In  case  2),  a  data  input  routine  (DATAIN)  is  executed.  Data  is  read  from 
the  current  data  set  and  compared  with  the  list  of  desired  inputs,  as  speci¬ 
fied  by  the  original  parameter  file.  If  a  match  is  found,  the  data  is  stored 
in  an  array  for  further  processing;  otherwise  it  is  ignored.  In  either  case, 
data  input  continues  until  the  required  data  values  have  been  read.  In  the 
event  that  the  start  of  a  new  data  set  is  encountered  before  input  of  the 
current  set  is  completed,  an  abnormal  exit  from  the  data  input  routine  occurs 
and  an  error  message  is  printed. 

Regardless  of  the  mode  of  exit,  the  search  routine  is  invoked  to  locate 
the  next  data  set.  As  before,  the  header  is  tested  and  any  of  the  three 
conditions  previously  mentioned  will  halt  the  search.  This  time,  however, 
case  2)  causes  a  repeat  of  the  data  input  routine  and  case  3)  causes  a  data 
processing  function  (MANIP)  to  be  performed.  Case  1)  still  produces  an  error 
message,  but  continues  to  the  data  processing  step  rather  than  terminating  the 
program. 

The  data  input  routine  also  tests  for  the  following  error  conditions  to 
the  FLAG  subroutine: 

1 )  Broken  sensor. 

2)  Overloaded  sensor. 

3)  Value  exceeds  upper  set  point. 

4)  Value  below  lower  set  point. 

All  of  these  conditions  are  indicated  by  flags  which  are  present  in  the 
original  data  from  the  data  logger.  These  tests  may  be  tailored  to  other  data 
formats  by  alteration  of  the  FLAG  subroutine. 

The  data  storage  array  is  organized  as  a  two  dimensional  matrix  in  which 
the  columns  contain  data  obtained  from  a  particular  channel  and  the  rows 
correspond  to  different  data  sets  (different  times).  MANIP  accesses  a  cross 
reference  matrix  and  determines,  for  each  input  channel,  the  column  in  which 
the  corresponding  data  has  been  stored.  This  column  is  processed  in  accor¬ 
dance  with  the  function  specified  for  that  input  channel  and  the  resulting 
value  ic  replaced  in  the  data  array. 

In  general,  the  function  will  be  different  for  each  type  of  sensor  and 
will  have  been  chosen  so  as  to  convert  the  data  logger  output  (typically  a 
voltage)  into  a  value  of  an  appropriate  physical  quantity  having  the  desired 
units.  During  this  processing  step  the  units  of  the  input  quantity  (as  read 
from  the  original  data  file)  are  compared  with  the  expected  input  units  (as 
given  in  the  parameter  file)  and  an  error  message  is  produced  if  a  disagree¬ 
ment  is  found.  This  message  identifies  the  date,  time,  and  channel  for  which 
the  error  was  detected  and  also  shows  what  units  were  actually  found. 
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After  processing  of  the  data  is  completed,  the  output  subroutine  (DATOUT) 
is  called. 

Table  3  shows  a  sample  of  the  input  data  obtained  from  the  data  logger. 
The  resulting  error  and  warning  messages  appear  in  an  output  file, 

FOR06.DAT.  To  this,  DATOUT  appends  a  tabular  listing  of  data  from  the  select¬ 
ed  input  channels,  as  shown  in  Table  4.  The  year,  Julian  date,  and  time  are 
listed  in  the  left  columns.  For  each  selected  channel,  a  column  of  data  will 
be  produced  having  a  heading  which  gives  the  channel  number,  sensor  identifi¬ 
cation,  and  the  units.  A  i,«ximum  of  twelve  channels  of  data  may  appear  across 
a  line  printer  page.  If  more  channels  were  requested,  additional  pages  will 
be  produced,  each  having  the  date  and  time  on  the  left  and  appropriate  column 
headings  across  the  top. 

In  addition  to  selection  of  the  data  to  be  listed,  the  user  formats  the 
output  by  providing  FORTRAN-type  format  specifications,  as  desired.  A  total 
of  ten  characters  (including  spaces)  should  be  specified  for  each  channel 
which  is  to  be  printed. 

An  additional  file,  FOR07.DAT,  is  also  produced  by  METSRT.  This  file, 
which  contains  all  of  the  processed  data  plus  reference  information  needed  by 
subsequent  programs,  was  intended  to  be  read  only  by  computer.  The  format  was 
chosen  for  compactness  and  few  concessions  to  human  readability  have  been 
made.  A  sample  of  this  output  is  shown  in  Table  5. 

III.  METCDC  OPERATION 

METCLC,  listed  in  Appendix  B,  reads  both  the  control  file  and  file 
FOR07.DAT  and  calculates  values  for  the  following  derived  quantities: 

1 )  Absolute  wind  velocity,  expressed  as  a  wind  speed 
(knots)  and  bearing  (degrees  referenced  to  true  north). 

2)  Atmospheric  moisture  loading,  with  water  vapor 
concentration  in  ppm  by  volume. 

For  each  calculation,  the  name  and  channel  number  for  each  input  and 
output  channel  is  printed  in  the  summary  listing,  FOR06.DAT  (Table  6).  Any 
problems  encountered  (missing  channels,  for  example)  are  also  listed  at  this 
point. 


The  calculated  values  are  then  stored  in  the  data  array  and  are  available 
for  output  in  both  a  tabular  form  and  as  an  array  intended  to  be  read  by 
subsequent  programs.  As  before,  any  combination  of  these  results  may  be 
selected  for  listing  by  setting  the  appropriate  print  switches.  The  print 
formats  are  specified  by  the  user.  An  example  of  this  listing  is  given  in 
Table  7. 

To  provide  versatility  and  make  future  alterations  simple,  calculations 
have  been  implemented  in  separate  subroutines. 

In  general,  there  may  be  several  alternate  sources  of  data  for  these 
subroutines  (multiple  anemometers  or  hygrometers  may  have  been  used,  for 
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Table  3 

An  excerpt  from  the  USNS  HAYR3  cruise  data  file  (FOR05.DAT). 
The  entire  file  is  over  200  pages  long  when  printed. 
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Table  4 

The  METSRT  listing  (FOR06.DAT)  generated  using  the 
parameters  and  data  filp  shown  in  Tables  2  and  3. 


No  warnings  or  error  messages  were  produced  by  this  data 
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Table  5 

Part  of  the  METSRT  data  dump  (FOR07.DAT)  produced  from  the 
inputs  shown  above.  This  data  serves  as  the  input  to  METCLC. 
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The  METCLC 


data 

run 


Table  7 

listing  (FOR06.DAT)  from  the 
which  produced  Table  6. 


urns  mi  cpoxse 

DUH  fOB  THE  l*il U  29  ,  198 


same 


TsaonGH  THO 


29, 


1981 


CHANNEL 

II 

rel.hun. 
?2?  CENT 

1981 


0 

77. 

’ 00 

77. 

200 

18. 

3  00 

77. 

i*00 

70. 

500 

78. 

600 

79. 

700 

30. 

800 

80. 

900 

80. 

1  000 

07. 

1100 

74. 

1200 

7a. 

1300 

73. 

*400 

7  3. 

*500 

73. 

1600 

73. 

1700 

73. 

1900 

72. 

1900 

72. 

3000 

6  9. 

2100 

69. 

22U0 

6  8. 

2300 

71. 

CHANNEL  CHBBB8X. 
20  100 

PRESSURE  ABS.5PD. 

TORS  knots 


758.  a 
759.9 

759.  0 
759.2 

759. “ 
759.1* 
759.5 
759.5 
759.5 

759.9 
760.1* 

760.  # 
760.  » 

760.  « 

761.  2 

761.9 
761.  5 
761.7 
761.7 

761.7 

761.8 

761.3 

761.4 
761. u 


10. 

13. 

9. 

5. 
8. 

6. 

7. 

6. 

8, 

11. 

9. 
15. 

3. 

15. 

19. 

19. 

12. 

12. 

15. 

11. 

10. 
8. 
8. 
9. 


CHUNKS!. 

101 

ABS.Oia. 

EEGRE2S 


89. 

81. 

80. 

82. 

82. 

88. 

93. 

99. 

91. 

89. 

95. 
81. 
60. 
62. 

7  0. 
51. 

96. 
59. 
9  9. 
41. 
26. 

7. 
349. 
34  0. 


uirtNN 8L 
102 

120  HAP. 
HU 


1  5426. 

1 i)*2. 

1  544  2. 
147a  1. 
14404. 

1  4261. 

1 11>  7. 

1 i 340. 
13211. 
tijJ2. 
11340. 
14612. 
14154. 
11111. 
14020. 

1  1140. 
11112. 
1j7»4. 
Ila27. 

1  4366. 
12/12. 
12762. 
12518. 
121/6. 


CHANNEL 

900 

BCXS  2 
PPSV 


.  i  6E*05 
.  1 6E  *05 
. 168*05 
. 16B*05 
.7  55*05 
.15S«05 
. 1 5E*05 
.  1  95*05 
. 145*05 
.  1 4E *05 
.  1  45  *05 
.155*05 
.155*05 
.1  55*05 
.155*05 
.155*05 
.1  5**05 
.1*5*05 
.1 »E*05 
.1  55*05 
.145*05 
.1  *C»05 
. 1 4E*05 
.13r*05 


2  FUNCTION  CALLS 


! 


12 


'W 


example).  Accordingly,  there  is  provision  for  user  specification  of  the 
inputs  for  each  calculation.  In  fact,  the  same  calculations  cam  be  repeated 
with  different  combinations  of  inputs  and  the  results  may  be  listed  for  com¬ 
parison. 

Absolute  wind  velocity  is  calculated  by  vector  addition  of  the  absolute 
velocity  of  the  sensor  platform  (ship)  and  the  relative  wind  velocity.  Sub¬ 
routine  input  and  output  vectors  are  in  the  form  of  a  magnitude  and  a 
direction.  It  is  assumed  that  the  direction  is  in  degrees  from  true  north 
(for  absolute  bearings)  or  degrees  clockwise  from  the  platform  velocity  vector 
(for  relative  bearings)  and  that  speeds  are  in  knots. 

Atmospheric  moisture  loading  is  calculated1  as 

P  (T  )  „ 

[H20]  =  H  104 

a 

where  [H^O]  =  water  vapor  concentration  (ppmv);  H  =  relative  humidity  (%); 

Pa  =  ambient  pressure  (mb);  Pg(Ta)  =  saturation  vapor  pressure  (mb)  at  ambient 
temperature  Ta  (°K). 

The  saturation  vapor  pressure  may  be  obtained  from 

P  (T  )  =  P-  exp  ["  y'* C  tnl 
s  a  0  n  -I 

where  PQ  =  1013.25  mb;  C,  =  13.3185;  C2  =  -1.9760;  C3  =  -0.6445;  C4  *  -0.1299 
and  t  is  given  by 


a 

with  T.  =  ambient  temperature  (°K). 

Inputs  to  this  subroutine  are  assumed  to  be  in  units  of  torr,  per  cent, 
and  degrees  centigrade  for  pressure,  relative  humidity,  and  temperature, 
respectively.  The  output  is  the  water  vapor  concentration  in  parts  per 
million  by  volume. 

IV.  METPLT  OPERATION 

METPLT  is  designed  to  plot  selected  subsets  of  the  data  contained  in 
FOR (3 7 •  DAT  according  to  the  specifications  given  in  file  F0RjJ4.DAT.  The  actual 
details  of  generating  a  plot  file  are  handled  by  DISSPLA,  a  package  of 
FORTRAN- callable  subroutines  provided  by  Integrated  Software  Systems 
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Corporation.  Basically,  METPLT  provides  the  data  needed  by  the  DISSPLA  rou¬ 
tines.  Appendix  C  contains  the  source  code  for  METPLT,  but  not  for  any  of  the 
DISSPLA  software. 

The  parameter  file  (Table  2)  is  searched  until  the  "PLOT  PARAMETERS" 
section  is  found  and  the  title  to  be  used  on  the  output  is  read.  Next,  the 
data  file  is  read  and  stored  in  memory.  In  the  event  that  either  of  these 
files  is  missing,  or  if  the  "PLOT"  section  is  not  found,  an  error  message  will 
be  written  into  FOR06.DAT. 

The  next  parameters  to  be  read  specify  the  number  of  days  for  which  the 
data  is  to  be  plotted  on  a  single  page  ( NDAYS )  and  the  dates  of  the  first  and 
last  data  which  is  to  be  plotted.  The  set  of  all  plotted  data  for  an  NDAYS- 
long  period  is  referred  to  as  a  plot  set.  Typically,  the  length  of  a  plot  set 
is  seven  days  so  that  the  plotted  output  will  have  one  week  of  data  per 
page.  There  may  be  more  than  NDAYS  between  the  initial  and  final  dates  speci¬ 
fied,  in  which  case  multiple  plot  sets  will  be  produced.  Each  plot  set  may 
itself  involve  several  pages  of  output  since  there  are  normally  only  three 
graphs  per  page. 

We  must  still  specify  which  data  is  to  be  plotted  and  how  it  is  to  be 
plotted.  This  is  done  by  providing  sets  of  parameters  which  define  each  axis 
for  each  plot.  Subroutines  XAXIN  and  YAXIN  are  responsible  for  reading  and 
storing  these  parameters. 

METPLT  first  searches  for  the  set  of  parameters  which  describes  the 
desired  X-axis,  then  it  looks  for  corresponding  Y-axis  parameters.  For  each 
X-axis,  there  may  be  multiple  Y-axis  specifications  so  that  several  different 
graphs  may  easily  be  generated  using  the  same  independent  variable. 

Each  axis  specification  consists  of  the  following  nine  parameters: 

1 )  Channel  number. 

2)  Channel  name. 

3)  Channel  units. 

4)  Minimum  value. 

5)  Incremental  value. 

6)  Maximum  value. 

7)  Threshold  value. 

8)  Hysteresis  parameter. 

9)  Axis  type. 

The  channel  number  tells  the  program  which  data  is  to  be  plotted  on  the 
specified  axis.  In  the  event  that  a  channel  number  of  zero  is  given,  METPLT 
will  use  time,  rather  than  data  values,  for  that  axis.  In  this  case,  the  axis 
will  be  labeled  with  the  Julian  dates  and  each  day  will  be  labeled  at  1200  and 
2400  hours.  For  purposes  of  axis  specification,  however,  times  must  be  given 
in  minutes. 

The  channel  name  and  units  are  used  to  produce  a  label  for  the  axis. 
Minimum,  maximum,  and  incremental  values  are  needed  in  order  to  calculate  the 
scale. 


Threshold  and  hysteresis  parameters  provide  increased  control  over  the 
plotted  output.  If  the  value  of  any  coordinate  is  below  the  corresponding 
threshold,  plotting  of  the  point  will  be  suppressed.  The  hysteresis  param¬ 
eters  allow  points  to  be  suppressed  if  they  lie  within  a  specified  “dead  band" 
surrounding  the  most  recently  plotted  point.  Note  that,  when  hysteresis  is 
set  to  zero  on  any  axis,  the  dead  band  area  will  also  be  zero  and  all  points 
will  be  plotted.  In  order  to  prevent  this,  any  of  the  hysteresis  parameters 
may  be  set  to  a  negative  value  and  will  then  be  ignored. 

The  axis  type  parameter  allows  the  user  to  select  either  a  linear  Y-axis 
(type  =  LIN)  or  a  “vector"  Y-axis.  Only  linear  x-axes  are  permitted  in  the 
current  version  of  METPLT. 

In  the  vector  mode,  the  two  channels  representing  R-  and  ©-components  are 
designated  by  type  =  RVEC  and  type  =  AVEC,  respectively.  They  are  used  to 
generate  a  vector  quantity  which  is  then  displayed  as  an  arrow  of  the  appro¬ 
priate  length  and  direction.  The  tip  of  the  arrowhead  is  located  at  the 
corresponding  X-coordinate  for  the  quantity.  For  reference,  a  short  arrow 
(not  of  unit  length)  is  drawn  in  the  zero  degree  direction  and  a  scale  is 
provided  on  the  Y-axis.  This  scale,  the  axis  name,  and  the  axis  units  are 
those  given  for  the  vector  magnitude  channel. 

The  vector  plotting  mode  is  useful  for  representing  quantities  such  as 
wind  direction. 

Provisions  have  been  incorporated  into  METPLT  allowing  easy  program 
enhancement  to  include  other  types  of  axes,  such  as  logarithmic  scales  or 
vectors  expressed  as  X-  and  Y-components. 

Error  messages  are  written  into  FOR06.DAT  by  XAXIN  or  YAXIN  if  the  axis 
type  is  not  defined  or  (in  the  case  of  vector  axes)  if  the  two  types  are  not 
self  consistent. 

At  this  point,  subroutine  SETUP  determines  the  first  and  last  dates  for 
the  next  plot  set  and  searches  the  ITIME  array  to  locate  the  corresponding 
rows  in  the  data  matrix.  If  they  cannot  be  found,  an  error  message  is  pro¬ 
duced  and  the  plot  set  is  skipped.  Assuming  that  the  rows  have  been  located, 
LOAD  copies  the  appropriate  data  into  the  X-  and  Y-arrays  needed  by  the  actual 
curve  drawing  routines.  During  the  loading  process,  each  datum  is  checked  to 
see  if  it  is  invalid  (the  character  string  " - "),  if  it  is  below  the  thres¬ 

hold,  or  if  the  hysteresis  criterion  is  not  met.  In  the  first  two  cases,  the 
point  will  not  be  plotted  and  a  message  to  this  effect  will  appear  in  the  plot 
summary  listing  (FOR06.DAT).  In  the  third  case,  a  message  will  also  be 
printed  but  the  point  will  not  be  suppressed  unless  the  hysteresis  test  fails 
for  all  axes. 

After  all  of  the  graphs  on  a  page  have  been  completed,  a  page  caption  is 
added.  The  title  (specified  in  the  parameter  file)  is  written  across  the  top 
and  a  subtitle,  giving  the  initial  and  final  dates  of  the  plot  set,  appears 
below  the  title. 
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Examples  of  the  printed  and  plotted  output  from  METPLT  are  shown  in 
Tables  8  and  9. 

Additional  pages  are  produced  as  necessary  in  order  to  graph  all  of  the 
data  in  the  first  plot  set.  Further  plot  sets  will  then  be  created,  each  one 
starting  on  the  Julian  date  following  the  end  of  the  previous  set. 

V.  METINP  OPERATION 

METINP  is  the  only  interactive  program  in  the  METEOR  package.  A  listing 
of  the  FORTRAN  program  is  given  in  Appendix  D  and  an  example  of  a  terminal 
session  is  shown  in  Appendix  E.  METINP  uses  this  dialog  to  construct  the 
parameter  file,  FOR04.DAT. 

Initially,  the  user  is  asked  to  identify  the  program  for  which  a  param¬ 
eter  file  is  desired.  The  answer  to  this  question  is  used  to  select  one  of 
three  major  subroutines:  SORTIN,  CALCIN,  or  PLOTIN.  These  produce  control 
files  for  METSRT,  METCLC,  and  METPLT,  respectively. 

In  the  case  of  CALCIN,  additional  information  is  requested  regarding  the 
specific  type  of  calculation  required.  Depending  on  the  response,  either 
WNDIN  (for  wind  velocity  calculations)  or  MSTIN  (for  moisture  loading)  will  be 
called. 

In  order  to  minimize  the  size  of  the  program,  each  parameter  input  rou¬ 
tine  utilizes  the  same  set  of  I/O  subroutines.  Subroutine  TTYIN  writes  a 
prompter  message,  reads  the  user's  response,  and  stores  the  answer.  Subrou¬ 
tine  FILE  then  formats  the  answer  as  required  for  that  specific  line  of  the 
parameter  file. 

To  make  these  two  subroutines  more  generally  useful,  they  operate  on  data 
arrays.  For  example,  the  prompter  character  string  and  the  corresponding 
input  format  are  contained  in  arrays  PROMPT  and  FORMIN,  respectively,  while 
the  user  response  is  stored  in  INARA Y.  Each  call  to  TTYIN  or  to  FILE  may 
therefore  be  tailored  to  specific  requirements  by  passing  the  appropriate 
arrays  as  arguments  in  the  subroutine  call. 

After  all  required  information  has  been  obtained,  subroutine  CHECK  writes 
the  parameters  to  the  TTY  in  exactly  the  form  in  which  they  will  appear  in  the 
final  parameter  file.  If  any  changes  are  required,  subroutine  EDIT  allows  the 
old  line  to  be  overwritten  by  a  new  one,  which  is  then  displayed,  when  all 
lines  have  been  verified,  the  complete  set  is  written  onto  the  disk.  If 
further  input  is  desired,  the  entire  process  may  be  repaated,  either  for 
another  function  for  the  same  program  or  for  a  different  program. 

Since  METSRT  is  the  first  program  to  be  used  in  data  analysis,  it  is 
assumed  that  a  new  file  FOR04.DAT  will  be  required  whenever  a  METSRT  control 
file  is  to  be  created.  For  this  reason,  SORTIN  causes  a  new  disk  file  to  be 
opened  and  any  existing  file  with  the  name  FOR04.DAT  will  be  destroyed.  Some 
care  must  therefore  be  exercised  to  ensure  that  the  current  file  (if  one 
exists)  has  been  saved  under  another  name  before  a  new  file  is  opened. 
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Table  8 

The  METPLT  summary  (FOR06.DAT)  resulting  from  the 
parameters  and  data  of  Tables  2  and  5,  respectively 
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Table  9 

METCLC  plotted  output  corresponding  to  the  above  summary. 
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METCLC  and  METPLT  both  require  the  output  of  METSRT,  so  it  is  assumed, 
whenever  control  files  for  these  programs  are  requested,  that  the  METSRT 
parameter  file  already  exists.  Accordingly,  the  new  parameters  are  appended 
to  the  existing  file,  which  is  not  lost. 

VI.  SUMMARY 

METEOR  provides  a  coordinated  set  of  programs  which  can  read  data  tapes, 
locate  specified  types  of  data,  test  for  a  wide  range  of  error  conditions, 
calculate  values  of  several  derived  quantities,  and  produce  both  printed  and 
plotted  output,  all  under  control  of  a  user-created  command  file.  Existing 
functions  may  be  selected  as  required  and  new  functions  may  be  added  with 
relative  ease. 

Although  originally  intended  to  process  meteorological  data,  this  soft¬ 
ware  package  should  be  equally  applicable  to  any  situation  in  which  large 
quantities  of  diverse  data  are  acquired  over  long  time  periods. 

In  many  cases, data  may  be  collected  on  several  different  data  logger 
systems  simultaneously.  For  these  situations,  it  would  be  advantageous  to  be 
able  to  merge  the  resulting  data  files.  Other  possible  improvements  include 
addition  of  statistical,  curve  smoothing,  and  cross-correlation  capabilities 
in  METCLC  and  provision  for  better  control  over  plot  size  and  shape  in  METPLT. 

We  expect  that  other  users  may  find  it  necesssary  to  revise  the  METEOR 
programs  to  meet  their  special  requirements.  It  is  hoped  that  the  documenta¬ 
tion  provided  will  prove  to  be  sufficient  for  this  purpose.  Any  comments  or 
suggestions  regarding  alterations  to  or  extension  of  these  programs  will  be 
welcomed. 
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Appendix  A 

Listing  of  Program  METSRT 
(Version  2.0) 
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Appendix  B 
*9  of  Program 
(Version  i.q) 
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Appendix  C 

Listing  of  Program  METPLT 
(Version  3.2) 
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Appendix  D 

Listing  of  Program  METINP 
{Version  1.2) 
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Appendix  E 

METINP  Terminal  Session 
Creation  of  a  Parameter  File 
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Appendix  F 

METSRT  Error  Messages 

Several  references  to  the  error  detection  capabilities  of  METSRT  have 
previously  been  made.  In  this  section,  these  capabilities  will  be  discussed 
in  more  detail  and  some  examples  will  be  give. 

Five  different  classes  of  errors  may  be  distinguished,  as  follows: 

1 )  Failure  to  find  some  piece  of  required  information. 

2)  Errors  occurring  during  data  input. 

3)  Conditions  which  were  detected  and  flagged  by  the 
data  logger. 

4)  Inability  to  properly  process  some  type  of  data. 

5)  Errors  arising  from  the  processing  of  a  specific 
piece  of  data. 

Errors  of  the  first  type  are  detected  in  the  main  program  when  an  attempt 
to  OPEN  a  file  fails  or  when  the  keyword  "SORT"  cannot  be  found  in  the  param¬ 
eter  file.  In  all  cases,  execution  is  immediately  aborted,  and  the  reason  is 
printed  out. 

Subroutines  SEARCH  and  DATAIN  look  for  errors  of  the  second  type,  which 
include  such  things  as  read  errors,  duplicate  channel  numbers,  and  incorrectly 
formatted  data  records.  Each  such  problem  is  listed  and,  insofar  as  possible, 
analysis  of  the  remaining  data  continues. 

Once  a  record  has  been  read  and  found  to  be  valid,  it  is  checked  by  FLAG 
to  determine  if  any  of  the  data  logger  warning  flags  were  set.  Broken  or 
overloaded  sensors  are  detected  at  this  point  and  these  data  are  ignored.  If 
the  data  logger's  preset  upper  or  lower  limits  were  exceeded,  this  will  also 
be  be  detected  but,  in  these  cases,  the  data  is  not  rejected.  Finally,  if  an 
illegal  character  appears  in  one  of  the  positions  reserved  for  data  logger 
flags,  this  fact  will  be  reported. 

In  general,  each  type  of  data  will  require  a  special  function  subprogram 
to  perform  the  appropriate  calculations  and  conversion  of  units.  MANIP  uses  a 
computed  GO  TO  statement  to  select  the  correct  function  for  each  channel.  In 
the  event  that  the  channel  number  lies  outside  the  range  of  the  GO  TO,  or  if 
no  function  has  been  defined  at  the  specified  statement  label,  then  an  error 
message  will  be  printed.  The  data  will  be  left  in  its  original  (input)  form. 

During  processing  of  a  specific  data  record,  various  errors  may  occur. 

For  the  most  part,  it  is  left  up  to  the  individual  function  subprograms  to 
make  any  tests  which  may  be  required.  For  example,  the  values  of  various 
power  supply  voltages  are  tested  and  messages  printed  if  they  lie  outside  the 
specified  ranges. 

One  general  requirement  of  all  data,  regardless  of  type,  is  that  it  must 
have  the  correct  units  (those  expected  by  the  corresponding  subprogram).  For 
each  channel,  the  proper  units  are  declared  in  the  parameter  file  at  input 
time.  DIDDLE  compares  the  actual  units  with  the  expected  units  and  lists  all 
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discrepancies . 

Most  of  the  error  messages  presently  incorporated  in  METSRT  are  illus¬ 
trated  in  the  following  three  tables.  Table  FI  shows  the  complete  error 
listing  as  it  appears  in  FOR06.DAT.  In  Table  F2,  those  messages  produced 
during  data  input  are  related  to  the  errors  in  the  data  file  which  caused 
them.  Table  F3  similarly  correlates  processing  time  error  messages  with  the 
erroneous  data  records. 

In  all  of  these  examples,  the  parameter  and  data  files  were  specifically 
constructed  to  exercise  the  maximum  possible  number  of  error  trapping  rou¬ 
tines.  Not  illustrated  are  those  (such  as  “NO  DATA  FILE  FOUND")  which  are 
mutually  exclusive. 
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Table  FI 

METSRT  Error  Messages 


iiUi  lESiAGt  DEHC. 

a*mci 


ISIS  307  2200  CH  A MEL  18  :  LOHEP  LIBIT  EX  CE  ID  ED. 

1979  307  2200  CHANNEL  22  :  SEISO?  'ifERLOAD 

1979  307  2200  CH  All  CL  23  :  UNEXPECTED  Ch  AS  ACT  EX  :  2  -i .  -  X 

1979  ?07  2200  CHAM  EL  2*  :  SEISO?  950KEI 

IICORPIBHEISI BLE  DATA  AT  2200  HOURS  Cl  DAI  307,197*; 

A  CAC  LIS l 

SISTEH  HESrt  AFTER  2200  BOORS  ON  DAT  307,  1979.  DA  IA  Jill.**  )MP»!D. 
ACQUISITECI  RCSOBBO  AT  2230  BOORS  Cl  DAT  307,1979 
ILLEGAL  RECORD  AT  2230  HOUSS  01  DAI  307,  1979  : 

////// 

1979  307  2230  CHAIIEL  17  :  UNEXPECTED  CHARACTER  :  *  " 

El  C  CP  LATA  S  FT  IOB  2230  HOORS  OR  DAT  307,1979  :  4  ..iAJSRLS  READ 

EIIOI  DO  SI  IG  SEARCH:  19  79  307  2230 
1001 979 

1979  307  2300  CHAIIEL  22  :  SEISOP  II  IAIGE 

1 S 79  307  2300  CHAIIEL  2«  i  SSBOI  IS  OK 

1979  307  2300  CHMSEL  23  *.  399EX  LIBIT  EXCEEDED. 

DOPLICATt  DATA  FCOID  F01  CHAIIEL  2»  AT  2300  HOURS  jj  3J 7,1 9' 9: 

A  29  -C.3100  t 

ECF  COIIIG  SEASCB:  1979  307  2300 

OISPFCIFIED  PROCESS  FOB  CHAIIEL  16  (ITERATION  ■  1) 

1  979  307  2200  CHAIIEL  18  :  HOHIDITT  OUT  OF  RANGE  l  -*.j;  CORRECTED  TO  0 

1979  3Q1  2200  CbAIIEL  21  :  HOHIDITT  OUT  OF  RAIGE  l  1  J  J.  |  ;  CORRECTED  TC  ’00 

1979  307  2200  CHAIIEL  23  :  OIEXPTCTPD  OUTS  (HI)  fJJ MJ 

1979  307  2230  CHAIIEL  23  :  CORRECT  OUTS  {  V)  F03  i  J 

1979  307  2230  CHAIIEL  25  :  H.  0.  TCLTAGE  OUT  OF  SAN*..  i  i.oi  7J  T| 

1979  JO’  2230  CHAIIEL  26  :  PCS.  ?IR.  SUPPLE  OUT  Of  iAfc-3  (  44.190  V) 

1579  307  2300  CH  ARIEL  27  :  REG.  I*  B.  SUPPLY  OUT  OF  (  0.030  »J 

ILLEGAL  CHAIIEL  IUHBEP  55  (ITEPATICI  •  13) 


Table  F3 
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Appendix  G 
Non-Standard  FORTRAN 

An  attempt  has  been  made  to  restrict  the  statements  used  in  these  pro¬ 
grams  to  the  set  defined  by  1966  ANSI  standard  FORTRAN.  However,  in  several 
places  non-standard  statements  have  been  used.  Often  this  was  done  because 
the  desired  function  was  sufficiently  complex  that  no  simple  alternative  was 
available.  In  other  cases,  the  non-standard  statements  were  considered  to 
involve  relatively  trivial  functions  which  could  easily  be  deleted  by  other 
users  without  detriment  to  the  overall  program  function. 

In  this  appendix  we  briefly  discuss  these  non-standard  features  and 
suggest  possible  alternatives  for  some  of  them. 

I.  PROGRAM  'name' 

This  statement  assigns  a  name  to  the  main  program  just  as  FUNCTION  or 
SUBROUTINE  are  used  to  designate  subprograms.  It  may  be  omitted  without 
affecting  any  program  functions. 

II.  PARAMETER  M=n 

PARAMETER  M  =  n  assigns,  at  compile  time,  the  value  'n'  to  the  constant 
'M'.  In  the  METEOR  package,  PARAMETER  statements  are  used  to  set  MDATA  (the 
maximum  allowable  number  of  data  sets)  and  MCHAN  (maximum  possible  number  of 
data  channels).  These  two  constants  are  then  used  to  dimension  many  of  the 
arrays  in  both  the  main  programs  and  in  the  subprograms.  A  PARAMETER  state¬ 
ment  must  appear  in  each  subprogram  in  which  MDATA  or  MCHAN  are  to  be  used. 

If  the  PARAMETER  statements  are  omitted,  then  each  occurrence  of  MDATA 
and  MCHAN  must  be  replaced  by  explicit  values. 

III.  OPEN/CLOSE 

These  statements  control  the  characteristics  of  the  files  used  for  input 
and  output.  The  following  arguments  may  be  used  with  OPEN  or  CLOSE 
statements : 


1 ) 

UNIT 

=  n 

Defines  the 

number. 

logical  unit 

2) 

DEVICE 

=  ' DSK' 

Specifies  that  the  device  is  a 
disk. 

3) 

ACCESS 

=  ' SEQOUT' 

Initializes 

device  for  write. 

=  'SEQIN'  Sets  device  for  read. 

=  'APPEND'  Sets  device  for  write  but  does 

not  initialize.  New  data  will 
be  added  to  the  end  of  the 
existing  file. 

4)  DISPOSE  =  'DELETE'  Delete  file  after  it  is  closed. 

=  'SAVE'  Save  file  after  close.  This  is 

the  default. 

5)  FILE  =  'filename'  Allows  new  files  to  be  named. 
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The  default  name  is  FOR0n.DAT, 
where  'n'  is  the  logical  unit 
number. 

6)  ERR  =  s  Causes  a  branch  to  statement 

number  's'  if  an  I/O  error 
occurs. 

In  many  systems,  the  functions  of  the  OPEN  and  CLOSE  statements  may  be 

performed  by  job  control  commands  external  to  the  program.  However,  the  error 

recovery  function  may  not  be  available  in  these  cases. 

METINP  closes  and  reopens  file  FOR04.DAT  at  several  points.  These  state¬ 
ments  could  be  eliminated  and  the  file  allowed  to  remain  open  continuously 
during  program  execution. 

In  METCLC,  subroutine  DATOUT  closes  FOR07.DAT  for  input,  reopens  it  for 
output,  and  rewrites  the  entire  file  (with  modifications).  The  equivalent 

effect  might  be  achieved  by  defining  a  new  logical  unit  to  receive  this  output 

[change  the  WRITE  (7,f)  statements  to  WRITE  (8,f),  for  example],  deleting  the 
old  file  7,  and  renaming  the  new  file  (file  8,  in  this  example). 

IV.  STOP  'string' 

This  statement  causes  the  message  'string'  to  be  written  to  the  default 
device  (TTY  for  interactive  jobs,  LOG  file  for  batch  jobs)  at  the  time  that 
the  STOP  is  encountered.  These  statements  serve  little  purpose  in  batch  jobs 
(in  most  cases  the  same  message  is  available  in  FOR06.DAT)  but  have  proven  to 
be  convenient  in  debugging  from  a  terminal.  They  may  be  replaced  by  standard 
STOPS. 

V.  RETURN  n 

This  statement  allows  subroutines  to  return  to  any  point  in  the  calling 
program.  Any  subroutine  that  uses  this  feature  must  have  one  or  more  ’&s' 
arguments  (where  's'  is  a  statement  number)  in  the  CALL  and  corresponding 
dummy  arguments  in  the  SUBROUTINE  statement.  A  RETURN  n  will  then 

return  to  the  statement  number  represented  by  the  nth  asterisk  (counting  from 
the  left). 

A  substitute  for  this  function  might  involve  setting  the  value  of  a 
variable  within  the  subroutine  and  then  using  a  computed  GO  TO  in  the  calling 
program  to  branch  to  the  desired  statement  number. 

VI.  END  =  s  /  ERR  —  s 

The  END  =  s  feature  is  used  as  part  of  the  READ  statement  to  direct  the 
program  to  statement  's'  if  an  End  of  File  (EOF)  is  read.  The  format  of  the 
statement  is  [READ  ( n , f , END  =  s)  'list'].  Since  this  is  only  used  to  enable 
the  program  to  print  out  an  appropriate  message  before  termination  of  the 
program  it  is  not  really  necessary  and  may  be  omitted. 

If  it  is  desirable  to  retain  this  feature,  the  function  could  possibly  be 


simulated  by  placing  some  standard  character  in  the  last  record  of  each 
file.  The  input  may  then  be  tested  for  the  presence  of  this  character  and  an 
EOF  routine  called  if  it  is  found.  This  procedure  might  be  implemented  by 
doing  a  READ  with  an  A-format,  then  using  a  BACKSPACE  and  another  READ  (with  a 
different  format)  if  the  special  character  was  not  found. 

ERR  =  s  is  also  used  in  a  READ  statement  in  essentially  the  same  fashion 
as  the  END  function.  ERR  operates  as  described  in  the  discussion  of  OPEN  and 
CLOSE  statements.  Since  the  error,  if  present,  is  detected  by  the  operating 
system's  I/O  routines  there  is  little  that  can  be  done  to  mimic  this  function 
within  the  FORTRAN  program.  However,  it  is  possible  that  the  job  control 
language  may  provide  commands  by  means  of  which  an  error  recovery  may  be 
accomplished. 

VII.  SKIP  RECORD  n 

The  SKIP  RECORD  statement  causes  the  next  record  on  device  ’n'  to  be 
skipped  during  input.  It  is  equivalent  to  the  construction: 

READ  (n, f ) 
f  FORMAT  (/) 

VIII.  ENCODE/DECODE 

These  two  statements  allow  data  to  be  reformatted  within  the  computer. 
They  both  require  arguments  as  follows: 

ENCODE  (n,f, 'array' )  'list' 

DECODE  (n,f, 'array' )  'list' 

where  'n'  is  the  number  of  characters  to  be  transferred  and  ’ f '  is  a 
format  statement  number. 

ENCODE  is  somewhat  like  WRITE  in  that  information  in  the  variables 
specified  by  'list'  is  transferred  to  a  string  under  control  of  a  FORMAT 
statement.  However,  instead  of  being  written  to  an  output  device,  the  string 
is  written  into  variable  'array'. 

Conversely,  DECODE  reads  'n'  characters  contained  in  'array',  formats 
them  as  specified  by  FORMAT  statement  'f,  and  stores  the  results  in  the 
variables  given  in  'list'. 

ENCODE  is  used  by  METPLT  to  create  character  strings  which  are  used  as 
captions  and  axis  labels.  These  strings  could  be  explicitly  defined  in  the 
program,  read  in  from  a  file,  or  they  may  be  omitted  entirely  without  signifi¬ 
cantly  altering  the  functions  of  these  programs. 

METINP  reads  a  record  once  in  an  A-format,  then  uses  DECODE  to  reformat 
the  string  as  required.  This  could  be  accomplished  by  using  a  BACKSPACE 
followed  by  another  READ. 


